arXiv:1502.04892vl [cond-mat.quant-gas] 17 Feb 2015 


Discrete Breather and Soliton-Mode Collective Excitations in 
Bose-Einstein Condensates in a Deep Optical Lattice with 
Tunable Three-body Interactions 

Galal Ahmed Alakhaly and Bishwajyoti Dey 
Department of Physics, University of Pune, 

Pune-411007, India. 


Abstract: We have studied the dynamic evolution of the collective exci¬ 
tations in Bose-Einstein condensates in a deep optical lattice with tunable 
three-body interactions. Their dynamics is governed by a high order discrete 
nonlinear Schrodinger equation (DNLSE). The dynamical phase diagram of 
the system is obtained using the variational method. The dynamical evo¬ 
lution shows very interesting features. The discrete breather phase totally 
disappears in the regime where the three-body interaction completely dom¬ 
inates over the two-body interaction. The soliton phase in this particular 
regime exists only when the soliton line approaches the critical line in the 
phase diagram. When weak two-body interactions are reintroduced into this 
regime, the discrete breather solutions reappear, but occupies a very small 
domain in the phase space. Likewise, in this regime, the soliton as well as 
the discrete breather phases completely disappear if the signs of the two-and 
three-body interactions are opposite. We have analysed the causes of this 
unusual dynamical evolution of the collective excitations of the Bose-Einstein 
condensate with tunable interactions. We have also performed direct numer¬ 
ical simulations of the governing DNLS equation to show the existence of the 
discrete soliton solution as predicted by the variational calculations, and also 
to check the long term stability of the soliton solution. 

1. Introduction 

The unprecedented control over ultra-cold atoms in optical lattices allows 
the Bose-Einstein condensates (BECs) formed by such atoms to model vari¬ 
ous physical phenomena. A few examples are Bloch oscillation, superfluid to 
Mott-insulator transition, Ehmov physics etc. Ehmov state were hrst exper¬ 
imentally observed in BECs [1] 35 years after its prediction by Ehmov [2]. In 
the Ehmov state of the condensate, the three-body interaction can completely 
dominate over the two-body interaction. For a very recent paper on the Eh¬ 
mov ehect in BECs see [3]. Ehective multi-body interactions has also been 
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experimentally observed in BECs |1]. It is well known that the Feshbach res¬ 
onance mechanism allows hne tuning of the two-body interactions in BECs. 
In a recent paper Safavi-Naini et. ah |5] proposed an interesting mechanism 
for tuning the onsite three-body interaction using Ehmov states EE] and its 
radio frequency (RE) coupling with the triply occupied state. They found 
that a strong three-body attraction affects the order of the Mott-insulator- 
superfluid transition. Such RE coupling of the trimer state has also been 
experimentally observed in ultra cold Fermi gases [7]. Very recently, Daley 
and Simon suggested an experimental method to make coherent three-body 
interactions dominate the physics of ultracold lattice gases [8]. 

BEC exhibit a range of many-body phenomena, such as collective excita¬ 
tions which are signihcantly affected by the inter-atomic interactions. Soliton 
and breather modes are examples of such collective excitations. These non¬ 
linear modes are localised in space, and arises due to the balance between 
nonlinearity and dispersion in the system [9] . Such localised collective modes 
has been experimentally observed in BECs [10]. These nonlinear collective 
modes also occurs in nonlinear lattice systems and are known as discrete 
breathers and discrete solitons. Nonlinear lattices are used to model energy 
localisation and transport in physics and biology [n]. Collective excitations 
modes of a fermionic gas of ®Li atoms in the BEC-BCS crossover regime has 
also been experimentally observed [12] and it shows a critical transition with 
variation of the strength of inter-atomic interactions. Analytical studies of 
the effects of two- and three-body inteactions, and also the anharmonicity 
of the trap potential on the stability and collective excitations of BECs, are 
also reported in the literature [13]. However not much is known about the 
localised collective excitation of BECs loaded in an optical lattice. It is nat¬ 
ural to ask how the two- and three-body interactions affects the nonlinear 
localised collective excitations of BECs in a deep optical lattice. 

In this paper, we study how the dynamical evolution of the collective 
excitations in BEC is affected by the presence of local and tunable three-body 
interactions. For this we obtain the dynamical phase diagram of the system 
using the time-dependent variational method. We hnd that the dynamical 
evolution of the collective excitations of the condensate in the regime where 
three-body interaction dominates over the two-body interaction shows very 
unusual and interesting behaviour. We have also performed direct numerical 
simulations of the governing equation of the system to show the existence 
of the discrete soliton solutions as predicted by the variational calculations, 
and also to check the long term stability of such solutions. 
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The plan of the paper is as follows: In Section 2, we derive the higher 
order discrete nonlinear Schrodinger eqnation (DNLS), which describes the 
dynamics of the scalar BECs in a deep optical lattice. In Section 3, we de¬ 
scribe the variational dynamics to stndy the phase diagram of the system. 
In Section 4, we report varions phases of the system as obtained from the 
variational analysis and plot the phase diagrams. In Section 5, we report re- 
snlts of the direct nnmerical simnlations of the higher order DNLS eqnation. 
Finally we give a conclnsion in Section 6. 

2. Model 

The dynamics of the BEC with two- and three-body interactions is given 
by the effective mean-held ID quintic time dependent Gross-Pitaevskii (GP) 
eqnation or the nonlinear Schrodinger eqnation (NLS) 

+ C/0 I $ P +C/2 I (1) 

where Qq and g 2 are the conpling constants for the two- and three-body 
interactions which depend on the s-wave scattering length a and the atomic 
mass. The condensate is trapped in an optical lattice potential given by 
Vext = Ui{x,y) sin^(27r2:/A) where A is the wave-length of the laser, A/2 is the 
lattice spacing of the optical lattice and Ui{x, y) is determined by the intensity 
of the laser beam. Eq. (1) has been derived from varions approaches snch as 
renormalisation gronp approach [T3], from many-body Schrodinger eqnation 
[15] and also from microscopic theory [16]. For the deep optical lattice (for 
large enongh laser power) one can nse the tight-binding approximation, where 
the strongly localised wave fnnctions (atomic orbitals) at the well of the 
lattice can be approximated by the Wannier fnnctions. Accordingly, we write 
the condensate order parameter $(r, t) as a linear combination of the atomic 
orbitals (LGAO) 0(r — r„) = 0ri(i') localised in the trap n as 

$(r,f) = •\/7^^V’n(t)0(r-r„) (2) 

n 

where Nt is the total nnmber of atoms. The orthonormal conditions of the 
atomic orbitals (Wannier fnnctions) gives 

J dr0„(r)0„+i(r) = 0 

/ dr<Pl{r) = 1 (3) 


3 


where is the n-th amplitude which can be written in terms of the 

number of particles Nn and phases 9n in trap n as 


■Ut) = ('i) 

where pn = Nn/Nx- For the tight-binding approximation (Eq. (2)) to be 
valid, the depth of the optical lattice should be much larger than the chemical 
potential p and the energy of the system should be conhned to the lowest 
band ira. Substituting Eq. (2) in equation Eq. (1), the GP equation reduces 
to the discrete nonlinear Schrodinger equation (DNLSE) with higher order 
nonlinearity [IHlCniED] 

1 

(V'n-l + ^n+l) + (^n + Aq | V'n P +A 2 | HV’n (5) 

at 2 

where 

Ao = J dr[0^], (6) 

(T) 

and 

^ ^ ^ ^ext4>l]- ( 8 ) 

Here time has been rescaled as t ^ where 

f 

A" J -1“ 0n^ea;t^n±l] 5 (9) 


Eq. (5) can be obtained from the Hamiltonian equation of motion 


V’n 


dH 


( 10 ) 


where the Hamiltonian function is given by 

H = + C^n+l) + | T + y I V’n Ay | V’n |®] (H) 


with *-0* and ijjn as conjugate variables. Both the Hamiltonian and the norm 
Sn I '0n P= 1 are conserved. 
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3. Variational Dynamics 


We study the dynamics of the BEC using the time-dependent variational 
method. The basic idea is to take a trial function which depends on a num¬ 
ber of time-dependent parameters, and to derive equations of motion for 
these parameters using the variational principle. The amplitude of the wave- 
function determines the density distribution, while its phase determines the 
velocity field. As mentioned above, the interaction between the particles 
leads to collective behaviour of the many particle system. It is therefore use¬ 
ful to write the variational wavefunction in terms of variational parameters, 
such as a center-of-mass coordinate and the width of the condensate, hence 
describing the collective behaviour of the condensate [2T] . For this, we study 
the dynamical evolution of a Gaussian prohle wave packet [T 8 l | 22 l |23] 

+ ip{n - C) + - C)^| (12) 


where the variational parameters C,{t) and 7 (f) are the center and the width 
respectively of the density pn =\ tjjn P, vif) and 5{t) are their associated 
momenta, and fc is a normalisation factor. The linear site dependence of 
the phase term shows that the velocity or the conjugate momentum of the 
center-of-mass coordinate (^{t) is site independent. Similarly the quadratic 
site dependence of the phase term shows that the velocity or the conjugate 
momentum of the width 7 (f) depends linearly on site n [21]. The dynamical 
evolution of the variational wavepacket can be obtained by a variational on 
the Lagrangian 

n=oo 

L= Y. (13) 

n=—oo 

Using Eqs. (11) and (12), and after rearrangement, L can be written as 


7^5 

Ao 


2 A 2 

8 

_20r7_ 


_3\/37r7^_ 


-|-cospe ^ —U( 7 ,C) (14) 


where 


and 


V 


1 7 ^ 5 ^ 


V = k J dn[ene 


(15) 

(16) 


5 










The equation of motion for the variational parameters g*, where g* = P, C, 'y 
and 6, can be obtained from the Euler-Lagrange equations 


as 


A 

dt \dqi 

P = - 


dqi 


dV 


dC 

( = sm{p)e~^ 

k / ^ x2\ -r? I 2 Ao I6A2 

t) = (— - d ) cospe ’ H-+ 


AdV 


T 




and 


7 = 7 ^ cos pe ^ 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


The pairs (C, p) and (^, 6) are canonically conjugate dynamical variables 
with respect to the effective Hamiltonian 


H = _ —h „ — cospe *^ + 1 ^( 7 , C) 


2\/^7 3\/37r7^ 

The wave packet group velocity is given by 

dH 
dp 


( 22 ) 


Vg = 


(23) 


and the inverse effective mass is given by 

1 d^H 


m 


dp 


— = cospe 


-V 


(24) 


The numerical solutions of the above coupled nonlinear ODEs (Eqs. (18- 
21)) along with the Hamiltonian (Eq. (22)) gives the dynamical regimes or 
phases of the system (Eq. (5)). Knowing the solutions, various phases can 
be obtained as follows: as mentioned above, 7 and <5 are conjugate variables. 
7 represent the width of the wavepacket and 6 denote its conjugate mo¬ 
menta. A plot of S(t) versus ■yit) gives the corresponding phase space (5— 7 ) 
trajectory. To obtain the soliton phase, we note that the solitons are shape¬ 
preserving nonlinear localised solutions. This means that the width of the 
soliton solution do not change with time. A soliton solution can be obtained 
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by imposing the condition 7 = 5 = 0. This implies that the soliton mode is be 
represented by the hxed point of the 5 — 7 trajectory. The discrete breathers 
are nonlinear collective excitations which are localised in space and oscillates 
in time. This implies that for the discrete breather solutions, the width 7 (t) 
and its conjugate momenta 6{t) will oscillate in time and the corresponding 
5 — 7 trajectories will be closed curves. Similarly, for the self-trapped phase 
the motion of the center of the wave packet stops 0 ) as the effective 

mass tends (Eq. (24)) to inhnity. On the other hand, for the diffusive phase 
there is complete spreading of the wave packet and this phase corresponds 
to 7 cx), C 7 ^ 0, therefore the effective mass is hnite. Various phases as ob¬ 
tained from the solutions of Eqs. (18-21) can be shown in the phase diagram. 


4. Phase diagram 


For a horizontal optical potential, the onsite energy e„, as well as V((C, a), 
are constant [18]. Accordingly from Eq. (18) the momentum p{t) is con¬ 
served. This implies p{t) = Pq where Pq is the initial momentum. From Eq. 
( 22 ) the conserved initial energy can be written as 


Ho = 


Aq 2A2 

—H-=-cos Poe 

2y/na 3^/3na 


8 


(25) 


For the sake of convenience we have denoted the width of the wave packet 
by a instead of 7 , where 0 = 7 ^. Eq. (25) can be rewritten as 


52 


8 q! log 




—Ho 


cos Po 


+ 4 




(26) 


The trajectories in the {a — 5) plane can be obtained from the above equation 
along with the solutions of the coupled equations, Eqs. (18-21). In the 
region cospo > 0) there are two phases. This can be seen as follows: since 
the argument of the logarithm term in Eq. (26) has to be positive, this 
implies that for Hq > 0 the width a{t) can attain a maximum value amax 
for f —)■ 00 . The value of amax can be obtained by equating the argument 
of the logarithm term to zero. When a —)■ a^nax, we obtain from Eq. (26) 
that 5 —)■ 00 . Eq. (15) then implies p —)■ cx 3 and accordingly we obtain 
from Eq. (19) that C ^ 0, and from Eq. (24), that 1/m* —)■ 0. This 
corresponds to the self-trapped phase since the centre of the wave packet 
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stops as the effective mass goes to infinity. Thus we see that even if the 
system is fully Hamiltonian (conservative), the apparent damping occurs 
due to the fast (exponential) growth of the effective mass with time. On the 
other hand, when — cospo < Hq < 0, we see from Eq. (26) that in this case 
a{t —)■ cxd) —)■ cxD, and accordingly 5 —)■ 0 and from Eq. (15) ?7 ^ 0. But 
from Eq. (19) C 7 ^ 0 and from Eqs. (24) and (25) we see that the effective 
mass is hnite as ^ ~ —Hq > 0. There is complete spreading of the wave 
packet, giving rise to the diffusive phase. The critical line which separates the 
self-trapped and the diffusive phases is obtained from the condition Hq = 0. 
Substituting TTo = 0 in Eq. (25), we get the critical line (with (5o = 0) as 

Ao,c = 201^ [cospoe~^/^^° ~ o 1 

\ dVoTrao/ 

The critical line for cospo > 0 is shown in Fig. 1. 

The negative effective mass (cospo < 0) case is more interesting. This 
is because in this case, beside the self-trapping and the diffusion regimes, 
we get two new phases. These are the discrete soliton and discrete breather 
phases. In this case, the conserved initial energy Hq and the equation for the 
(a — (5) trajectory are same as Eq. (25) and Eq. (26) respectively, except 
that the cospo terms in these two equations should be replaced by — | cospo|. 
Following the same procedure as discussed above, we can see that in this 
case there are also self-trapping and diffusion regimes, and the critical line 
separating these two regimes is obtained from the condition IIq = | cospo|. 
Replacing Hq by | cospol in the expression of the conserved Hamiltonian, we 
obtain the critical line which separates the two regimes as 

Ao,c = cos Pol (1 - ) (28) 

\ 3v dTrcto / 

As mentioned above, the soliton mode is represented by the hxed point of 

the a — 6 trajectory ( d = 0 and 5 = 0). Eqs. (20) and (21) then have 

stationary solution a(t) = ag and 6(t) = 0 (a(t) = j^(t)). From Eqs. (15) 

_ 

and (19), we obtain that ( = sin(po)e = constant. This implies that for 
soliton solutions, the center of mass moves with constant velocity, and the 
width (shape) of the soliton do not change with time. From Eq. (20) we 
then obtain 

Ao,soliton = 2 ./^ (I cospo|e-'/ 2 “° - 

V tto V SvStt / 
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(29) 




Figure 1: Dynamical phase diagram with two- and three-body interactions 
(A 2 = 1) with ao = 2, (5o = 0 and Co = 0. 



This gives the soliton line in the phase diagram. The soliton line gives the 
value of the two-body interaction parameter Aq for which the soliton solu¬ 
tions are allowed. It may be noted that the soliton line corresponds to the 
maximum of the energy (Eq. (22)). This is because Eq. (29) can also be 
obtained from the condition ) = 0 and for which ) < 0. 

The soliton line is shown by the dotted line in Fig. 1. Yet another inter¬ 
esting phase which appears in the case of the negative effective mass is the 
discrete breather. Discrete breathers are nonlinear collective modes which 
are localised in space and periodic in time. For discrete breather solution 
a{t) , 5{t) and m* oscillates with time. The center of mass travels with 
nearly constant velocity and with an oscillating width. Discrete breather 
solution appears when the nonlinearity parameter Ao,c < Aq < Aq, soliton- The 
trajectories in the a — 6 plane are closed and a oscillates between the initial 
values tto and the maximum value > aQ. The complete phase diagram 
for the case when the two-body interactions dominates over the three-body 
interactions is shown in Fig. [H 

We can now obtain the phase diagrams for the case when the three- 
body interaction dominates over the two-body interaction (A 2 >> Aq). For 
this, we first consider the case when the three-body interactions completely 
dominate over the two-body interactions, and accordingly we set the two- 
body nonlinear interaction parameter Aq = 0. By substituting Ao,c = 0 in 
Eq. (28), we get the critical line A 2 ,c which separates the self-trapped and 
the diffusion regimes as, 

A 2 c ^ 3\/37rao _ g-i/ 2 «o^ (' 3 Q) 

Similarly, by substituting Aq, soliton = 0 in Eq. (29) we get the soliton line as 

A 2 ,soliton = Icospol (31) 

The soliton solutions are expected when the three-body interaction parameter 
A 2 values are on the soliton line, i.e. A 2 = A 2 ,soliton- From the above equation 
we can see that this can be satisfied for arbitrary value of cto (initial width 
of the condensate density). We checked this prediction of the variational 
calculations by direct numerical simulations of the governing equation of the 
system, the DNLS equation (Eq. (5)). However, from the direct numerical 
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Figure 2: Dynamical phase diagram of the BEC with only three-body inter¬ 
action (A 2 > 0 and Ao = 0) with ao = 100, (5o = 0 and Co = 0. The solid line 
represent the critical line (Ao,c) and the overlapping dotted line represent the 
soliton line (Aosoiiton)- The inset show the phase diagram for cospo > 0. 
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Figure 3: a — 6 trajectories in the state with only three-body interaction 
(A 2 > 0 and Ao = 0). the hgures are plotted for different values of A 2 > 0 
and cospo < 0 with ao = 5, po = 37r/4. The critical values of the three- 
body interaction parameter for different transitions are A 2 ,c = 2.7462 and 
A 2 ,soiiton = 2.61112. Inset shows the a — 6 trajectories for A 2 very close to 

A2,soliton- 




integration (Section 5) of Eq. (5) we obtain an unusual result that the discrete 
soliton mode collective excitations exist only for the large value of Oq. From 
Eqs. (30) and (31) we see that for large value of cto, A 2 ,soiiton —t A 2 ,c, i-e. the 
soliton solutions exist only when the soliton line approaches the critical line in 
the phase diagram. Fig.[2]shows the phase diagram for large value of ao = 100 
and we can see that in this case the soliton line (dotted line) approaches the 
critical line (solid line). From the direct numerical solution of the DNLS 
equation we hnd that indeed, there is a stable soliton solution for ao = 100. 
Details are given in Section 5. Yet another unusual result for the BEC 
with only three-body interactions is that it does not allow discrete breather 
collective mode excitations. The absence of the discrete breather modes in 
this regime can be seen from the following: discrete breather collective modes 
occur when the three-body interaction parameter A 2 satishes the relation 
A 2 ,c < A 2 < A 2 ,soliton- As mentioned above, the trajectories in the a — 6 
plane are closed, and a{t) oscillates between the initial value ao and the 
maximum value a™^^ > ao. The breather region extends until the value 
of parameter A 2 approaches A 2 , breather- For A 2 , soliton < A 2 < A 2 , breather, ao 
oscillates between ao and the minimum value a™™ < ao- Both the values 
a™“ and a™™ are roots (together with ao ) of Eq. (26) for Ao = 0 and 
cospo < 0. The condition for which this equation does not have another 
root a™"" < ao gives Aq, breather [H]- After rearrangement, one can show that 
A 2 , breather E giveu by the maximum value of the r.h.s of the expression 


A2 

A2,soliton 


2aoa; 
(1 -x) 


_ g(a:-l)/2«oa: 


(32) 


Here a“™ = xao and x < 1. From the maximum value of the r.h.s, we find 
that A 2 ^breather t A 2 ^soliton- ^YheU A 2 ^breather t A 2 ^soliton, the area euclosed 
within the a — 6 trajectories shrinks to zero (to the soliton hxed point) 
implying nonexistence of the discrete breather modes. Also, for the three- 
body interaction parameter in the range A2,c < A 2 < A2,soliton there are no 
closed trajectories in the a—6 plane and this implies that there are no discrete 
breather solutions. This is shown in Fig. [3l where we can see that there are 
no closed trajectories around the soliton fixed point. The inset of this figure 
shows no closed trajectories very close to A 2 , soliton- However, if we now add 
a weak two-body interaction (Aq < A 2 ), and if both the parameters Aq and 
A 2 are positive, then the discrete breather solution is allowed. In this case, 
after rearrangement, we obtain 


13 






0.20 


0.15 


0.10 


0.05 

5 

0.00 


- 0.05 


- 0.10 


0 5 10 15 20 25 30 35 40 

a 


14 

Figure 4: a — 5 trajectories in the state with weak two-body interac¬ 
tion (A 2 > Aq) for different values of A 2 > 0 (beginning from the bottom 
A 2 = 0.5,1,1.5, A2,c, 2.32, 2.34, 2.36, A 2 ,soiiton, 2.382, 2.384, A 2 ,breather, 4,10) for 
the case cospo < 0 with cto = 5, po = 37r/4,Ao = 0.1. Inset 
shows the a — 5 trajectories close to A 2 ,breather (from the top A 2 ,c = 
2.39, 2.386, A 2 breather, 2.384, 2.383, 2.382). 








( 33 ) 


A2 

^2soliton 


2anx 
(1-*) . 


^{x-l)/2aox 


x/^x/x(l-x/x) . 

y/w{l—x) cospoe~^^^‘^0 0 


1 j_A 

2-^/7? cos poe“^/^“0 0 


Again, from the maximum value of the r.h.s of Eq. (33) we hud that 
A2,breather > A 2 ,soiiton, thus Satisfying the condition for the appearance of 
the discrete breather solutions. This is shown in Fig. |U where we can see 
that there are closed trajectories (on the right of the soliton hxed point) rep¬ 
resenting the discrete breather solutions. For A 2 > A 2 ,soiiton, but very close 
to A2,soliton, the area of the trajectories near the hxed point is very small but 
nonzero (shown in the inset of Fig. H]), thus allowing the discrete breather 
solutions to exist. Fig. O shows the complete phase diagram. 

When Ao have an opposite sign to that of A 2 , then both the solutions, 
soliton as well as the discrete breather, are not allowed. The absence of the 
soliton and breather modes in this case is due to the fact that A 2 ,soliton is 
much smaller than that of A 2 ,c and accordingly the A 2 ,soliton (soliton line) 
lies deep inside the diffusion region. Due to the exponential dependence of 
A 2 ,soliton on cto (Eq. (31)), it is not possible to increase cto till the A 2 ,soliton line 
cross the A 2 ,c critical line for the soliton solution to reappear. Similarly, the 
discrete breather solutions are also not allowed in this case since A 2 ,breather 
(breather line) also lies well below the A 2 ,c critical line. 


5. Numerical solutions of the higher order discrete nonlinear Schrodinger 
eqnation 


We have also solved the higher order DNLS equation (Eq. (5)) directly 
using numerical procedures. This verihes the existence of the discrete soli¬ 
ton solutions in presence of only three-body interactions, as predicted by the 
variational analysis. We have also checked the stability of the soliton solu¬ 
tions over a long time scale. We write tjjnif) = \J Pn(^)e*^"^*^ in terms of its 

two components, as 'ipn = {xn,yn) and 'ipn = {xn,yn)- Here and denote 
the real and imaginary parts respectively of In terms of these variables, 
the DNLS (Eq. (5)) can be written as the coupled equations 

Xn = —2 (l/n-l + yn+l) + ^ny-n + Ao(x^ -|- + ^ 2 {x‘^ -f- y^^yn (34) 

and 

yn '^(^Xfi—1 T Xn+l) ^nXn Ao(x^ -\- y^Xn A2(x^ -f- J/^) Xn (35) 
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Figure 5: Dynamical phase diagram in the state for positive two- and three- 
body interactions with ao = 2, (5o = 0, Co = 0 and weak two-body interaction 
parameter Aq = 0.24. Inset shows the small discrete breather region. 




Since we want to compare the numerical results with the variational results, 
we take the initial condition as the variational wavefunction (Eq. (12)) and 
accordingly write the initial values of (a;„, Hn) as 


X, 


= Vke 




cos 


p(n-C) + |(n-C)^ 


(36) 


and 


= Vke 


(r^-cr 


sm 


p{n-C) + ^{n-Cf 


(37) 


where a = V ■ We solve the coupled nonlinear equations (Eqs.(34, 35)) using 
the Runge-Kutta method, with the initial conditions given by Eqs.(36, 37). 
The constants of motion, the Hamiltonian (Eq. (11)) and the norm were 
checked at every time step of the numerical integration to verify that these 
remained constant over time. To connect the numerical results with that 
of the variational results, we calculate relevant parameters which shows the 
existence of the soliton solutions. Thus we calculate the numerical average 
position which is dehned as 


<n>= '^n\'tprV 

n 


(38) 


It is easy to show that < n >= ({t). As mentioned above (Eq. (12)), ({t) is 
the variational parameter which denote the average position or the center of 
the density =| p. We also calculate the numerical width of the wave 
packet 2 V<n^ > which is dehned as 

< V >= VrV — < n >‘^ (39) 

n 

Again, it can be shown that < >= 7^(f)/4 = Q((f)/4. As mentioned above, 

a{t) is variational parameter which denote the width of the wave packet (Eq. 
(12)). For the soliton solution to exist, the width < v? > should not change 
with time. Fig. [6] shows a stable soliton mode collective excitation as ob¬ 
tained from the direct numerical integration of the DNLS equation (Eq. (5)) 
for ao = 100. From the hgure we can see that the shape of the soliton solu¬ 
tion do not change with time. The inset shows the average position and the 
width over a long interval of time. This conhrms the long term stability of 
the discrete soliton solution. 
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6. Conclusion 


In conclusion, we have studied BECs in a deep optical lattice with a tun¬ 
able three-body interaction. We have obtained the dynamical phase diagram 
of the system analytically, using the time-dependent variational method. The 
phase diagrams shows interesting dynamical evolution of the collective exci¬ 
tations with variation of the strength as well as the sign (attractive or repul¬ 
sive nature) of the two- and three-body interactions. We have explained the 
causes of the appearance of the phase diagrams of the system. We have also 
solved the dynamical equation of the system, the higher order discrete non¬ 
linear Schrodinger equation (DNLSE), numerically, to check the existence 
of the soliton solution in the presence of only three-body interactions, as 
predicted by the variational analysis. As mentioned above, BECs where the 
three-body interactions dominate over the two-body interactions have been 
prepared experimentally. It should therefore be possible experimentally to 
look for the dynamical evolution of the collective modes of BECs as discussed 
above. 
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